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Abstract 



We compute structure constants in M = 4 SYM at one loop using Integrability. This 
requires having full control over the two loop eigenvectors of the dilatation operator for 
operators of arbitrary size. To achieve this, we develop an algebraic description called the 
O-morphism. In this approach we introduce impurities at each spin chain site, act with 
particular differential operators on the standard algebraic Bethe ansatz vectors and generate 
in this way higher loop eigenvectors. The final results for the structure constants take a 
surprisingly simple form. For some quantities we conjecture all loop generalizations. These 
are based on the tree level and one loop patterns together and also on some higher loop 
experiments involving simple operators. 
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1 Introduction 



One of the main reasons for trying to solve a gauge theory such as M = 4 SYM is to search for 
hidden structures pointing towards alternative descriptions of Quantum Field Theory. As 
a byproduct, we also hope to develop new techniques that will allow us to do computations 
which are currently beyond reach and which will improve our understanding of String Theory 
and Quantum Gravity. 

To illustrate how such structures are often hidden we will consider the problem of com- 
puting correlation functions of local gauge invariant operators in the planar limit of A^ = 4 
SYM. For simplicity, let us consider operators made out of two complex scalars only, say Z 
and X. As usual we identify them as spins for notational convenience, 

Tr{ZX...)^\U...} ■ (1) 

We use N for the number of X fields and L for the total number of fields. Then, the gauge 
invariant operators with definite total dimension A = L + 7 are the eigenvectors of the 
dilatation operator which can be written as a spin chain Hamiltonian as [1, 2] 

L 

H = J2 [( V ~ % 4 ) H n , n+1 + 2g 4 H n , n+2 + 0(g 6 )} . (2) 

n=l 

Here H nm = I nm — P nm , with I, P the identity and permutation operator, respectively. The 
energies of this Hamiltonian are precisely the anomalous part of the total dimension, 7. We 
can already identify remarkable simplifications which come about when considering physical 
quantities but which are very much obscured in the intermediate steps: take for example the 
dilatation operator (2). If we were to write its next few loop orders it would require pages 
and pages with no obvious structure present. Nevertheless, its eigenvalues, which are the 
physical quantities 7 can be obtained very neatly as we now briefly review. We can think 
of the states diagonalizing (2) as N excitations 4- moving with momenta p = {pi, . . . ,Pn} 
on top of a ferromagnetic sea of fields. The scattering of two magnons is governed by the 
magnon S-matrix S{pi,pj)} The quantization of these momenta is given by Bethe equations 2 

Lp k + y~] t logS" (puPk) = 2nn k , n k e Z . (3) 

l^k 

The S-matrix is known to all loop orders [4, 5] and is quite simple. Once the momenta are 
fixed the energy of the state is obtained by summing the energies e(pj) of each individual 
excitation, 

TV 

7 (p) = <Pj) > <P) = 8 9 2 sin 2 (p/2) - 32g i sin 4 (p/2) + 0(g 6 ) . (4) 

1 The scattering of N particles is also governed by the two magnon S-matrix since the many body S-matrix 
factorizes into a product of two body scattering processes in an Integrable theory. 

2 See section 3 for a very brief review and [3] for a more extended one. 
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This scattering picture needs to be corrected if the chains are too small [6, 7]. That is, the 
Bethe ansatz (3) is asymptotic. Equations (4) and (3) are a striking example of how beauty 
emerges at the end of the day for physical quantities. 

Once the dimensions of operators are known we know their two point functions. The 
next natural observables are the structure constants which govern the three point functions. 
For three point functions we also need the eigenvectors of the dilation operator. After all we 
need to Wick contract them (at leading order) and then decorate this contraction by loops. 
However, the eigenvectors can be horrendous at first sight. More specifically, at two loops 
(and higher) there will be some local corrections to the wave functions, called contact terms, 
when two or more magnons are close together. When four excitations are close together the 
corresponding two loop contact term correction is as messy as to occupy one full page!, see 
appendix A. 

For the spectrum problem we could avoid dealing with these contact terms by considering 
asymptotically large chains where there is always enough room for the excitations to move 
freely. Then, to quantize the momenta of the excitations all we care about is the S-matrix and 
to compute their energy we can consider a region where they are well separated. However, 
when we move to three point functions the contact terms come into play. When we contract 
the three operators we need to know the exact wave function. All regions will be probed 
and all of them will contribute. On the other hand, the structure constants C123 are physical 
quantities and as such, there ought to be some nice description where all the complications 
due to the contact terms and the loop corrections conspire together to give a neat final result. 

Confirming this expectation is the main goal of the current paper. To arrive at this 
result we propose an algebraic description of the eigenvectors of the dilation operator (2) 
which automatically incorporates all possible contact terms. This approach is based on the 
so called O-morphism which gives the name to the paper. 

Let us end this introduction by mentioning that there are numerous other exciting hints 
at hidden structures, other than the O-morphism [8, 9, 10, 11, 12]. They strongly inspired 
this work. It would be wonderful to explore in detail the different connections and differences 
between all these approaches. 

This paper is organized as follows. In section 2 we translate the problem of computing 
structure constants into the spin chain language. The rest of the paper is divided into two 
main parts. In section 4 we follow the coordinate Bethe ansatz approach combined with 
some inspired guesses based on our experience with the spectrum problem. This section will 
contain several conjectures for the one loop correction to the structure constants and also 
for some all loop generalizations. In sections 5 and 6 we follow an algebraic Bethe ansatz 
approach and we develop the O-morphism map. This will allow us to derive the most general 
SU (2) structure constants and prove all the one loop conjectures of section 4, based on the 
coordinate Bethe ansatz approach. We conclude in section 8. Several appendices contain 
complementary material to the main text. 



5 




-Eg 2 -^12+1 



O: 



Figure 1: Internal loop diagrams results in the insertions of the Hamiltonian densities at 6 
junction points [13, 14, 15]. This leads to (5) and (6) in the spin chain language. 

2 Structure Constants and Quantum Spin Chains 



In this paper we will consider structure constants in the SU(2) setup introduced in [3]. Let 
us recall how these objects can be efficiently computed using the spin chain picture. At tree 
level all we need to do is to Wick contract the three operators. These operators are the 
eigenvectors of the dilatation operator (2). At one loop level (a) we need to Wick contract 
the operators which are the eigenvectors of the two loop dilatation operator and (b) we need 
to include Hamiltonian insertions coming from Feynman diagrams [13, 14, 15] (see figure 1). 
In total 



where 



C 

norms 
simple 

involved 



123 



norms x simple x involved + 0(g 4 



(5) 



(6) 



V / J L 1 L 2 L 3 / V / (1|1)(2|2)(3|3), 
(i_^t ••• t\l-g 2 M N3 , N3+1 -g 2 U L3!l \3), 

N 3 L 3 -N 3 

(1|1 -g 2 U Ll2M2+l - g 2 m Lul \i 1 ...i Ll2 l^J) 

N 3 

(ii . . . i Ll2 t^t |1 - g 2 M Ll2 , Ll2+1 - <7 2 H L2jl |2) . 

L 3 -N 3 

Where i a can be ^ or J, and one sums over all possible 2 Ll2 intermediate states. Here |a) 
stands for the corrected two loop Bethe eigenstate corresponding to the operator O a . The 
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length and number of excitations of this operator are denoted by L a and iV a respectively. Of 
course, for these expressions to make sense and to be non zero we need 

N 1 = N 2 + N 3 and L 12 = L x - iV 3 = L 2 - L 3 + N 3 (7) 

where L 12 is the number of contractions between operators 0\ and 2 - For any pair of 
operators Oi and Oj we have 2Ljj = Li + Lj — Lk where k ^ i,j. 

The computation of structure constants in M = 4 SYM in the SU{2) sector up to one 
loop order is thus mapped into a very precise spin chain problem: we simply need to be able 
to compute all spin chain scalar products and expectation values appearing in (6). At order 
g° this was solved in [3, 16]. One of the main goals of the current paper is to compute the first 
quantum correction, of order g 2 . This would be the first step toward an all loop asymptotic 
expression for the structure constants. Obviously we want to find reasonable conjectures for 
all loop results based on the tree level and one loop results based on our experience with the 
spectrum problem. As an example, we will now consider the first contribution in (5), the 
norms. For that, we first need to work out the precise form of the two loop eigenvectors |a) 
which will allow us to ask nonambiguous questions about their norms. 



3 Two Loop States and Contact Terms 

The dilatation operator in M = 4 SYM can be thought of as a long-range spin chain Hamil- 
tonian whose interaction range coincides with the order of perturbation theory. At two loops 
we have (2). The single magnon state is given by 

T, "2-1 



\p) = e ipn \n) , |m, na, . . .> = 1 1 ■ ■ ■ t+t • • • tit • • •} = a™*™ • • • |t • • • t) • (8) 



n=1 ni-l 



This follows from translation invariance and is true to all orders in perturbation theory. The 
two magnon state is more interesting. At two loops we have 

|Pi>2>2>= J2 ^ Pini+iP2n2 + ^ P2ni+mn2 S(P2,Pi)){l + S^^ (9) 

ni<-«2 

The physical picture is clear: we have an incoming wave (e* Pl ™ 1+ip2ri2 ) describing two particles 
with momentum pi and p 2 . They can then collide. In 1 + 1 dimensions when two particles 
collide their momentum can at most be exchanged. This is a consequence of conservation of 
both energy and momentum. The outgoing plane wave (e ip2Ul+ipin2 ) describes the particles 
with their momenta exchanged. The relative coefficient between the incoming and outgoing 
plane waves is what we define as the S-matrix, which turns out to be given by [8, 6, 5] 

This is all we have as long as the particles are well separated. When they are within the 
interaction range of the Hamiltonian the two magnons will be in a complicated interacting 



7 



state. This is taken into account by introducing the so called fudge factors [6] which we also 
denote as contact terms. The purpose of this section is to study these terms. At two loops, 
and for two magnons, we only need 

C<3(P2,Pi) = -4sin(p 1 /2)sin(p 2 /2)sec((p 1 + p 2 )/2) . (11) 

The superscript indicates that this is a two loop effect while the subscript indicates that the 
two magnons need to be close to each other. 

The higher loop picture should be clear; new fudge factors have to be added with longer 
and longer range. For example, at three loops, the two magnon the state (9) is modified by 

the simple replacement of ^1 + 5 ni+ i >n2 g 

(i+«w.i [/ci 2 ;(p 2)Pl )+/ci 3 ;(p 2 ^ (12) 

Also, the S-matrix needs to be corrected further and at four loop order a new contribution 
to the S-matrix - the so called BES dressing phase [5] - shows up. For now let us continue 
our two loop discussion but move to a higher number of excitations. 

Interesting new structures arise when we move to three (and more) magnons. For three 
magnons we have 

\Pl,P2,P3) = 22 [ < ^123+02135 , 21+01325'32+03125 , 3l5'32+023l5 , 3l5 , 21+032l5'32S'3iS'2l]|ni,n2,n3) 
ni<«2<«3 

where S ab = S(p a ,p b ) and 

(p abc = exp(ip a ni + ip b n 2 + ip c n 3 )(l + g 2 5 abc ) . (13) 
Finally - and most relevant for the current discussion - 5 abc are the contact terms. We have 



'abc 



for ri2 — rii > 1 and n 3 — n 2 > 1 

C| 2 ,' (p a ,Pb) for n 2 - Hi = 1 and n 3 - n 2 > 1 

C^J(p b ,p c ) for n 2 — n\ > 1 and n 3 — n 2 = 1 

d?J,(p a ,p b ,p c ) for n 2 -n 1 = l and n 3 - n 2 = 1 



We see that on top of the two magnon contact term we now need a new contact term for 
the case when the three magnons are all close to each other. States with more magnons 
follow the same pattern. For magnons we will have A! terms in the wave function and 
A — 1 different contact terms: C„ (pi,p 2 ), C^J,(pi,p 2 ,p 3 ), . . . , Ci?).,(pi, . . . ,Pn)- By explicit 
diagonalization of the two loop dilatation operator we can compute the first few contact 
terms. We find 

C^(pi,p 2 ) = see equation (11) 

„(2i/ \ 2(cos(pi +p 2 ) + cos(pi +p 3 ) + cos(p 2 +P3) - 2cos(pi +p 2 +p 3 ) - 1) 

Q.. Pl,P2,P3 = : 1—, ^ ^ - f : : s 

1 - cos(pi) - cos(p 2 ) - cos(p 3 ) - cos(pi + p 2 + p 3 ) 

c i«.(Pi 5 P2,P3,P4) = see appendix A (15) 
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The complexity seems to grow in a scary and uncontrollable fashion. As we will derive later 
through indirect means, all the contact terms can be written very explicitly as ratios of 
complicated sums, see appendix C. In practice for N > 4 it takes forever to compute these 
sums analytically even using a computer (numerically one can easily go up to much larger 
iV's). For now, let us continue our experimental investigation and turn our attention to 
more physical quantities such as the structure constants, in the hope that the huge mess in 
the contact terms will somehow miraculously drop out. As argued in the introduction this 
should be the case and indeed it will be. 

Finally, let us end by recalling that the momenta of these Bethe eigenstates is quantized 
through the Bethe equations (3) which follows from the periodicity of the many particle wave 
function. The energy of the states is given by (4). 



4 Coordinate Bethe Ansatz. Educated Guesses. 

To compute one loop structure constants we need control over the spin chain scalar product 
and expectation values such as the ones in (6). We know them at tree level, that is for the 
one loop eigenvectors of the Dilatation operator. The purpose of this section is to guess their 
higher loop generalization by doing experiments. In particular, these experiments will allow 
us to find an analytic expression for the simplest nontrivial structure constant (corresponding 
to the case when 0\ is BPS). 

We denote the momenta of the Bethe roots of Oi, 2 and C 3 as k, p and q respectively. 
Sometimes we shall use |2) and sometimes we will instead write |p) etc. 



4.1 All Loop Norm Conjecture 



The simplest scalar products appearing in our list (6) are the norms of the Bethe eigenstates. 3 
At one loop the norm of a Bethe eigenstate is given by Gaudin's determinant [17, 3] 



(pIp) 



1 loop 



n 

A<3 



1/2 



det 

l<j,k<N 




^ + £-iog^°)( 



PhPk) 



0(A) (16) 



where p = {pj} are the momenta of the N magnons and L is the total length of the operator. 
At one loop the S-matrix is given by the leading order in (10) which we denoted by S^°\ 
Finally p* is the complex conjugate of pf, the first factor in (16) is absent for real momenta. 
Note that the expression inside the determinant is nothing but the left-hand side of Bethe 
equations (3) at leading one loop order. A natural guess for the higher loop generalization 

3 Usually a norm is not a physical object; it depends on the normalization of the state. In the case of the 
Bethe sates, however, the normalization is fixed uniquely by the clear physical picture of scattering of the 
magnons. Namely, we normalize the state by fixing to 1 the coefficient of the incident wave, see e.g. (9). 
After this the norm is well defined and is entitled with some physical meaning which is to be uncovered. 
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is that we simply have to correct the S-matrix 

1/2 



(pIp) 



all loops 



n 

i<3 



S(p*,P*) 



det 

l<j,k<N 




Lp k + log S{p h p k ) 



l^k 



(17) 



At first this guess may look oversimplified. After all, when going to higher loops we also 
needed to include the contact terms (15) which were quite complicated and one could suspect 
that the above result is correct modulo these contact terms. A simple way to settle this 
question is to perform some experiments. 4 We computed the norm for several different 
eigenstates with different number of magnons by brute force up to 4 loops. Remarkably, 
we find that (17) works perfectly when the contact terms are included! On the contrary, 
it is because of the very complicated contact terms that the final result turns out to be 
so remarkably simple! This clearly indicates that some deeper structure must be hidden; 
unveiling it (at least at two loops) will be the purpose of later sections. Remarkably, at 4 
loops the dressing phase is nontrivial and affects the eigenstate. Nevertheless we find that 
(17) works perfectly! This leads to the conjecture: 

For large enough operators, in the usual asymptotic sense, (17) holds to all loops. 



4.2 Double Vacuum Decay Amplitudes and Structure Constants 

We now have control over the norms contribution in (5), see previous section. We will now 
move to the remaining contributions: simple and involved. In this section we will 

• Compute simple 

• Compute involved for the simpler cases where 0\ is a BPS operator. 

It turns out that for 0\ BPS, the contribution involved is as simple as the simple term. 
More precisely, these quantities turn out to be related to an important scalar product which 
we call the Double Vacuum Decay (DVD) Amplitude, see below. Guessing the higher loop 
expression for this scalar product is the main goal of the current section. 



4.2.1 The Double Vacuum Decay (DVD) Amplitude 

We will now introduce an important building block of the structure constants which we call 
Double Vacuum Decay (DVD) amplitude. 

To define it let us describe zero energy states of the dilatation operator. We will work 
with a particular basis of zero energy states which are the so called vacuum descendents of 
length L and N spin flips. We denote them as 

/ T,\ 

permutations . (18) 




L-N 



Latter we will derive (17) for two loops from an algebraic Bethe ansatz approach where no guesswork 
will be involved. 
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We can glue two such vacuum states one after the other making what we call a Double 



Vacuum 



1 N 1 


1 N ' 1 






1 N ' 1 


1 L \ 


1 u \ 


>- 




1 L' 1 



The DVD amplitude is defined as a scalar 



product between such double vacuum and a Bethe state 

A L (q) = (-l) 



\N 



N'-N 1 


1 N 


L 


\ L'-L 



q 



(19) 



Here V is the length of the Bethe eigenstate |q) and N' is the number of magnons in it. If 
all momenta q are nonzero the Bethe eigenstate is a highest weight state. In this case the 
DVD amplitude does not depend on A. 5 

The DVD amplitude (19) is an important building block for the structure constants as 
we will now illustrate. For example, the simple contribution in (6) in the new notation 
becomes 



simple 



1 Ns 1 


I o 1 


1 N 3 \ 


1 La-Ns \ 



i-g 2 



W 3 ,N 3 +1 



+ e 



(20) 



Note that the bra state is annihilated by BL^+i unless + 1) = (N 3 ,N 3 + 1) or (L 3 , 1). 
This means that we can replace ii 13i _t 13+ i + Hl^i by X^i-^v+i m (20). In this way we 
recognize the one loop dilation operator, see (2), up to a factor of 1/2. Then we act with 
this operator to the right, on |3). At this loop order, we simply get (one half of) the energy 
of the Bethe eigenstate 



N :i 



7(q) = E 8 



S 2 sin 2 | + 0(^ 



(21) 



where q 



{qj} are the momenta of the magnons of the operator (9 3 . We conclude that 

1 



simple 



1 



7(q) Av 3 (q) 



(22) 



Physically, Al describes the probability amplitude for all the N magnons in |q) to occupy 
homogeneously the first L sites of the chain which has total length L' > L. Above N = N3 
and L = N3 as well which means so that magnons are all squeezed together in the first A^3 
sites. The scalar product AfXq) was studied in detail in [18] and will be considered again in 
section 4.2.3. 



4.2.2 Three Point Functions from DVD Amplitudes 

In the previous sections we computed the first contribution norms in (5) and reduced the 
second contribution, simple, to the computation of the DVD amplitude Ac,(q) defined in 
(19). We now turn our attention to the involved contribution in (6). 

In general, it is not possible to get rid of the Hamiltonian insertions in this quantity as 
easily as we did for the contribution simple. However, if 0\ is a BPS state then we can still 
do it. In this case, we can follow the steps of the previous section to obtain 

involved = M - i 7 (p)J ^-^(p) (23) 

5 The Bethe state is annihilated by S + . At the same time the adjoint operator S~ acting on the bra 
state would mix states with two various N from which we conclude that they must be equal up to a sign. 
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where p = {pi, . . . ,Pn 2 } are the momenta of the magnons of the operator (9 2 . For details 
see Appendix B. Combining this result with (22) we conclude that for 0\ BPS we have 

01 BPS: Cffifrq) = yfZj^ ( 1 - 2M±2M) A N M)A Ll . N3 ( P ) 

V 2 ; v /(^)( 2 |2)(3|3) 

O ia nd0 3 BPS: C5S(p,q) = y/Lj^U-^] ^^M = . (24) 

V 2 ; 4 /(£)(£a )(2 | 2> 



We already found the norms in section (4.1). We now turn to the DVD amplitude A. 



4.2.3 Conjectures for the DVD Amplitude 

In this section we find an explicit expression for the DVD Amplitude defined in (19) valid at 
tree level and at one loop. It is the last missing piece to obtain the full one loop corrected 
structure constants of the type (24). At tree level - i.e. for one loop Bethe eigenstates -we 
have [3] 6 

e n if" 

. , . aUa={p} fcea.fcSa feci /oe^ 
•At ree level (PJ = ^ • (25) 

3=1 l<i<j<N 

The sum is over all possible partitions of the set of momenta and |a| stands for the number 
of elements in the partition a. For example, for two particles, we have 

1 - f( PuP2 )e^ - f(p 2 ,Pi)e lL ^ + e iL ^ + P*) 

■Atiee level \Pl, P2) = T"— lV _ ivo TVTf s ■ 

(e 1 pi - l)(e IP2 - l)f(pi,P2) 

Finally, f(p, k) = 1 + i [u(p) — u(k)} _1 where u(p) was defined in (10). We will repeat the 
game of section 4.1 and try to guess what the higher loop generalization of (25) could be. 
The first natural guess is that this expectation value still takes the form (25) where we simply 
need to correct f(p,k). To check this proposal we computed (19) for a two-loop corrected 
operator O2 with two magnons. For that we use the two loop wave functions, with contact 
terms, as described in the previous section, see appendix D.l for details. We realize that 
(19) indeed still takes the form (26) with a corrected f(p, k) — > f(j>, k), 

ujp) - u(k) + % ( g 2 n( A C\ (07 , 

«"•*>= «&,)-«(*) I 1 ^(p)' + })(«(*)> + l) +ote) J- < 27 > 



(> 



As shown by Kostov, the sum (25) can be recast as a determinant [19], 
E II f(W)J[e«* = i< <teJu(p j ) k - 1 -e i * L (uM 

aUo=jp) kea,k£a feect 



fc-11 
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We should emphasize that this is far from obvious. In fact, it is because of the particular 
form of the contact terms that the result can still be put in the form (26). By computing 
(19) for a few magnons we verify that 

e (- i ) H n f(M)ji e<ts 

Mp) = *u«=M keaMa fea (2g) 

n(e-^-i) n f(R.Pi) 

j=l l<i<j<N 

works perfectly at one loop! In fact, by appropriately correcting (27) we seem to be able to 
reproduce the two and three loop results [20]. This leads to the next general conjecture: 

For large enough operators, in the usual asymptotic sense, (28) holds to all loops provided 
(27) is properly fixed. 7 

For more details on how to fix f(p, k) see appendix D.2. In particular, in this appendix, 
we present an alternative procedure which constrains this function to all loop orders. 

This concludes our study of the simplest correlation functions in (24) where 0\ is BPS. 
We conjectured the form of all the building blocks appearing in this expression and even 
proposed some all loop conjectures. 

On the other hand we should conclude this section by mentioning that we failed in 
guessing most general structure constants when both 0\ and O2 are non-BPS operators. 
This is partly because the tree level result simplifies enormously on-shell i.e. when the set 
of roots p is constrained by the Bethe ansatz equations. At the same time this complicates 
enormously the guesswork since this constraint is rather complicated to deal with when 
doing experiments. In sum, the guessing strategy used so far seems insufficient for these 
more complicated cases. We need some more scientific and systematic approach to find the 
one loop (and higher loops) deformation. Developing such an approach is the main goal 
of the current paper and will be the focus of the remaining sections. The more impatient 
readers can check Appendix F for the final result for the most general one loop structure 
constants for three non-BPS operators. 



5 Algebraic Bethe Ansatz Approach and Impurities 

In the previous section we followed the coordinate Bethe ansatz approach where the states 
are written as sums of plane waves with relative coefficients given by ^-matrix elements. 
Beyond one loop we also have contact terms, see section 3. In this section we will only 
consider the one loop states for which there is a well developed alternative approach known 

7 At higher loops the Dilatation operator has scheme dependent terms which cannot be fixed by symmetry 
and locality, see [21]. These terms have unfixed coefficients and do not affect the spectrum. They can be 
removed by a similarity transformation. The DVD amplitude is not invariant under such transformations and 
can, in principle, contain some scheme dependence which of course should cancel out in physical quantities. 
This subtlety will be considered elsewhere. 
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as the Algebraic Bethe ansatz. In this approach multi-particle states are simply given by 



l*;u>= E 

Non Intercepting Paths 
Moving Left and Up Only 



1 ttt 


mm 
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1 1 
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MAT 



W2 







L 



u-O-i/2 



(29) 

In fact, the states represented here already contain a new twist: the impurities 9j at each 
of the spin chain sites. To get the one loop Bethe eigenstates we set all impurities to zero, 
8j = and consider |0; u). However, as we will see in the next sections, these impurities will 
play a central role when going to higher loops. 

The relation between the Bethe roots Uj and the momenta pj is given by Uj = u(pj) where 
the change of variables was introduced in (10). This figure is a graphical implementation of 
the algebraic relation 8 

\6;u) = B( Ul )...B(u N )\t...t) • (31) 

It is quite remarkable that such creation operators B(uj) exist at all. Even though the 
different B's enter independently they create fully interacting particles! For example, for the 
two particle state, we obtain the following wave function (we have already set 9j = in this 



3 More explicitly, 



B(«) = <t|Q0 Ujo + 



u-0j- i/2 



io II) 



(30) 



Here the subscript jO indicates that the operator acts on the tensor product of physical Hilbert space Vj = C 2 
corresponding to spin chain site j and an auxiliary space Vq = C 2 . The bra and ket in this formula live in 
this auxiliary space so that B(u) is an operator acting on the full spin chain Hilbert space H = V\ <8> • • • ® Vl- 
Graphically, the action of many such operators can be represented as in (29). For algebraic Bethe ansatz 
reviews see for example [3] and [16]. 
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figure): 



111 + V 2 V f U 2+ 



U\ — i/2 J \u 2 — i/2 



^ ni m n 2 



11 2 



ttttlttttlttt 

L 




m=ni 



Y i f u 1+ ii2 \ m - ni f u 2 +ii2 \^- m - i f i Y"""Y ' V " 

Ui-i/2J u 2 -i/2 \m-il2) \u 2 -i/2) \u 2 + i/2J \u 2 -i/2) 

(32) 

which leads, up to a simple proportionality factor, to the exactly the same wave function 
i^coor{ni, n 2 ) found through the Coordinate Bethe Ansatz approach, 

Vw(n!,n 2 ) = e ^" 1+ ™ + S(p 2 , Pl )e ip ^ + ^ , (33) 

see (9) with g 2 = 0. Of course, this continues to hold for more magnons. For iV magnons, 
the relative normalization between the one loop coordinate states |p) on eioop and algebraic 
states |0; u) is given by 

N N 

|0;u) =^|p)onoioo P , H = Y[(e- ip i - l)Y[f{pj,Pk)- (34) 

j=l j<k 

where f(p, k) = l + i/[u(p) —u(k)]. Obviously, at the end of the day, all normalizations drop 
out for physical quantities. 9 



5.1 Impure States and Structure Constants 

To go to higher loops we need some sort of deformation of (31) consistent with Integrability. 
The simplest one, which basically costs us nothing, is to introduce so called impurities 6j at 
each lattice site j. In fact, in the previous section we already introduced these parameters, 
see (29). At one loop we simply set those impurities to zero. The idea - which will be 
explored in detail in the following sections - is to use these impurities to mimic the loop 
corrected propagation and interaction of the magnon excitations. This deformation will 
amount to acting on the states (29) with a differential operator made out of d/d9j and only 
then setting all impurities to zero at the end. 

Concerning the structure constants, this procedure will allow us to promote tree level 
results with impurities to higher loop results without impurities. For future use let us quote 

9 As we emphasized above the coordinate Bethe ansatz normalization is particular nice from a physical 
point of view and we would like to be able to derive at least the first several loop orders for the norm. In 
(17) we conjectured a value for the norm at all-loops in the coordinate Bethe ansatz normalization. 
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Figure 2: Tree level impure structure constant. Each box represents a sum of paths as in 
(29). We should also divide by the normalization of the three operators/spin chains. 

here what the three point functions at tree level look like once we introduce the impurities, 
see figure 2. This was worked out in [16], see also [22]. 

Let us make a simplifying assumption that O3 is BPS. In any event, the contribution of 
O3 factorizes into a simple contribution - see section 2 - so it is straightforward to generalize 
to the case with (9 3 non-BPS at any later stage if needed. 

We use Uj = u(pj) and Vj = u(kj) for the rapidities of operators 0\ and 2 respectively. 

Finally we have to introduce some notation for the impurities. We denote the impurities 
of each operator O a by 9j. Each impurity is associated to a line going from one operator 
to another so it is shared between two operators. We denote the set of impurities shared 
between O a and O b by 9 {ah) so that, for example, # (2) = 6 {12) U6 {23) and so on. When omitting 
the superscript it is implicit that we refer to the first operator, 6j = 9j. 

Then we have 

I ^tree with imp. i 

I 123 I / _ " / / /1 / -1 \ 1 /1 / 1 \ \ / /1 /q\ 1 /i/f)\ r" V / 



y/(0W;u\0W;u) (^);v|^);v)' 
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where 



In l/i \ T-f U k~ u m + i , , 9 lv^r. Ufc — O + i/2 1 -cA U k - U m - i 

{9; u\9; u) = det - — - > log — + - > log 

A i u k -u m j,k<NidUj i u k -9 a -%2 % *-f u k -u m + i 

(36) 

with a similar expression for the other norm. Before presenting the numerator in (35) let us 
discuss (36) briefly and make contact with the previous sections. The expression (36) is the 
norm of Bethe eigenstates with impurities. The Bethe equations with impurities take the 
form 

E l - log h + 'eg f^^i) = 2,„, , n.eZ (37) 

^ z \u k -6 a -i 2 *—f % \u k -u m + ij 

a=l x ' 7 m^k x 7 

which reduces to (3) if we set all impurities to zero. Hence the determinant in (36) resem- 
bles very much (16). There is however a different prefactor which arises because (36) was 
computed using the Algebraic Bethe Ansatz normalization whereas (16) was obtained in the 
coordinate Bethe Ansatz formalism. Of course, in (35), the normalizations of the kets are 
irrelevant. However, to compare the intermediate quantities (36) and (16) we better know 
the relation between the two. When there are no impurities the relative normalization is 
(34). Taking it into account we see that (36) and (16) agree in the homogeneous limit 9j — > 0. 



Finally we have the numerator which can be represented graphically as in figure 3. We 
have [16] 

N 3 JVi N 3 N 2 

n iik - ® + */ 2 )/ n - © + *m 

(0<V|<W> = ^ Wa det © 



n («™ - «») n ^ - v ^ n ^ - & 



n<m n<m n<m 



where 9$ = 9^ +1 _ m and 



(38) 



F nm 1 , Gnm \\ 



±» 9^+i/2U^ n (n k -v m +i) n^nK-^-i) 



V 



Where F nm is iVi x iV 3 and G nm is Ni x N 2 matrices completing together to a JVi x N\ square 
matrix. Note that only the 6» (1) = 9 {l ^ U 6>( 13 ) appears in the numerator factor. There is 
no dependence on the impurities 9^ associated to lines going from 2 to 0%. From the 
combinatorial point of view this is obvious, see caption of figure 3. 

The homogeneous limit is a bit more tricky here, see [16] for details. One has to resolve 
a 0/0 uncertainty using standard l'Hopital's rule. On the one hand the prefactor in (38) 
diverges as 9j — > 0. On the other hand the determinant vanishes since F nm will no longer 
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Figure 3: Main scalar product appearing in the impure structure constants (38). When 
contracting the three operators the rightmost sites of the spin chain O2 are contracted with 
empty sites in the chain of the third operator. Hence the weight associated with the lower 
right corner in the right picture is just 1 since the weight for each vertex where a blue thick 
lines crosses a thin gray line is 1, see (29). Because of this, this scalar product does not 
depend on the impurities 9^ at all and can be reduced to the object in the right [16]. 



depend on the column index m as 6 m — > 0. 10 

In what follows we will also need a particular case of the general scalar product (35). That 
is the case when 0\ degenerates to BPS. This can be done by sending all roots Uj — > 00. In 
this limit the state |1) becomes a vacuum descendent and (38) reduces to the generalization 
of the DVD amplitude with impurities, 



1 N!-N 3 1 


1 i 


1 Li-TVg 1 


1 L 3 -N 3 \ 



qUo=v 



L12 - 

a 



nn 



0„ + 'i/2 



a£an=l n ' J a£a ,a£a 



n /< 



a, a 



(40) 



e n - i/2 

where f(u,v) = l+i/(u — v). Denoting the square bracket by e^ 12 and using the same kind 



10 



Once we take the limit carefully the result is proportional to the determinant 



dT(v! 
dui 



dT(v t ) 
du N 



dT(v N2 



92 (ui) 

9T{vn 2 ) , , 



Qn 3 +i{ui) 



q n (u) = i{u + i/2)- n+1 - i(u + i/2)- n+1 



T{u) 
Q(u) 



1N1 



u-i/2 \ 1 Q(u+i) 



u+i/2 



Q(y-) 



(39) 



n,=i( M -%) 



The expression (39) is the most complicated factor in the tree-level structure constant for three non-BPS 
operators. In section (4) we tried to guess the higher loop generalization of several simpler scalar products 
rather successfully. Guessing the generalization of (39) turns out to be quite non-trivial as mentioned before. 
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of short-hand notation from [3] we can compress this expression to 



v) = ^ (-lpr^ , where (vac| = ( h^d ■ ( 41 ) 



(vac 

aUa 



Equation (35) - and also equations (36) and (41) - will be used in the next section. As we 
will see later, by use of well designed differential operators we will be able to promote tree 
level "impure" expressions to higher loops "pure" results. Developing the formalism to do so 
will be the subject of the next section. 



6 0— morphism and Quantum Corrected States 

In this section we introduce the main object of this paper: the 0-morphism. 

6.1 One Loop G— Morphim. Definition and Properties 

The O-morphism is a coupling dependent functional acting on functions of a set of variables 
6 = {9i}f =1 which we denote by double brackets 

6 : f{9) ^ (/(*)),- (42) 

In section 6.5 we will explain how one can define this operator by requiring it to obey 
several defining axioms. The justification of the name will be clear from what follows. At 
leading order the O-morphim is given by the differential operator 



(43) 



where 

V i = d i -d i+1 and di = ^-. (44) 

we use a periodic boundary conditions so that V L = &l — d\. For now let us take the 
definition (43) as our starting point. For further motivation see the next subsections, [22] 
and [20]. In particular, we will see below that by acting with this operator on the impure 
states of section 5.1, we generate perfect two loop eigenvectors of the Dilatation operator 
(2) including all the messy contact terms discussed in section 3 and appendix A! This is of 
course the main justification for (43). 

For now, let us explore some properties of the map (43) that will come in handy later on. 
Acting on a product of two functions we have 

f(6)h(6)) = (f(9)) (h(0)) + (f(8)]\h(9)) . (45) 
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where 

L 

=S 2 2V*. ( 46 ) 

In the right hand side we set to zero all impurities after taking the derivatives, just like we 
did in (43). In what follows this is always implicit; we always set 9j — > at the end. We 
denote the contribution (46) by cross-term. If it were not for this term then (45) would be 
the mathematical definition of a morphism 11 . Hence, the cross-term is the deviation from 
the morphism property. 

The cross-term (46) does vanish if one of the functions is a totally symmetric function of 
the 9i. That is we do have the restricted morphism relation 

Morphism : (f sym (e)h(e)) = (/ sym (0)) (h(9)) . (47) 

Another nice property that will be used shortly is that 

Dispertion : (V log = L log *£±$ + Otf) , (48) 

which one can easilly check directly using (43). In this expression x(u) is the so called 
Zhukowsky variable 

<u) + — = - , x(« = - - - + 0(g 3 ) . 
x(u) g g u 

These relations are very important and have interesting implications. For example, an obvi- 
ous consequence of the morphism relation (47) is that for any function F we have 

Um(0))) e = F[(f sym (e)) e ] . (49) 

Another nice implication follows when we combine both properties: Suppose we are dealing 
- as we often will - with a function f(u,6) which is also a function of the Bethe roots 
u = {u\, . . . , un}- These are quantized according to Bethe equations (37), 

V^l, (u k -6 a + i/2\ (u k -u m -i\ 
^ % \Uk-6a-i 2J t-f.i \u k -u m + ij 

a=l x ' ' myth v ' 

In particular, they depend on the impurities 9j in a totally symmetric way. Then, a conse- 
quence of the above relations is that | /(u, 9) J is still given by (43) where 



1. we can treat the Uj as independent of 8j when taking the derivatives in (43) but 



ii 



We are grateful to B. Vicedo for proposing this name to us 
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2. after applying the morphism - i.e. in the right hand side of (43) - the Bethe roots are 
now quantized according to the BDS equations [8] 

,u_ L , x{u k + i/2) Al fu k -u m -i\ 

<f>k = ~ log —, rj^r + > t log — = 2im k , n k e Z 51 

i x{u k -i2) i \u k -u m + ij 

instead of (50). The equations (51) are nothing but the corrected Bethe equations (3) 
and (10). 

Of course, such a prescription is exactly what one would like when dealing with expressions 
involving Bethe roots! 



6.2 Single Derivative of Monodromy Matrix et al 

The most important object in the Algebraic Bethe ansatz construction is the monodromy 
matrix 



L(9;u) = I 10 + — Pio ® • • - ® UiQ + ~ a ^Plo (52) 

Operators with subscript ij act on Vi ® Vj where Vi, . . . , Vl are the Hilbert spaces at each 
spin chain site while Vq is an auxiliary space which is also C 2 . That is, we can think of 
L(n) as a 2 x 2 matrix whose entries are operators acting on the full physical Hilbert space 
H = V 1 ®---®V L , 

mu) = ( #»> 6 («) ) (53) 
V C{u) D{u) J 

The operators B(u) are creation operators; the Bethe eigenstates are simply given by (31). 
Graphically, this action can be represented as in (29). A very important relation that we 
will use extensively is then 

L 

iV 3 L = L] + 5 jtL ^[H n ,„ +1) L] . (54) 

n=l 

where Hjj + i = Ijj+i — Pjj+i- As usual, we have set to zero all impurities after taking the 
derivatives. Hence, when acting on Bethe states \9; u) = B(u\) . . . B(un) It • • • t) with a 
single derivative Vj we get from (54) that Wj\9; u) = — Hjj + i|0; u) + S^l Z~2n=i HU,n+i|0; u) 
and since they are Bethe eigenstates the last term can be replaced by the state energy so 
that 

iVj\9] u) = -% j+1 |0; u) + 5 jtL T u \0; u) . (55) 

where 

N 



r " = E^r- ( 56 ) 



,=1 U i + 4 

Another very important object is the transfer matrix 

f(9;u) ee Tr L(9;u) = A{u) + D(u) . (57) 
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To the order we will be working at - i.e. as an expansion in small 9/s - this object commutes 
with the Hamiltonian Yln=i ^n,n+i- Hence, from (54) we find 

iDjf = -[% i+1 ,T]. (58) 

The identities (54), (55) and (58) will be extensively used in what follows. 



6.3 Integrability and Two Loop Eigenstates 

Consider the impure states (31). They diagonalize the transfer matrix (57), 12 

f(9;v)\9;u)=T(9;v)\9;u). (60) 

The eigenvalue T(9;v) is a symmetric function of the 9/s. A most striking fact about (60) 
is that the states \9; u) do not depend on the spectral parameter v and still they diagonalize 
the full family of operators T(9; v) for any value of v\ Integrability, that is exact solvability 
of the model, lies upon this remarkable fact. In this section we want to understand what 
happens to such remarkable properties once we apply the B-morphism to relations like (60). 

Let F(v) be an arbitrary function of T(9; v) and F(v) its eigenvalue. Then 

F(v)\9;u) =F(v)\9;u). (61) 

Applying the 0-morphism we get 

(P) (\9;u)) +(f]\\9;u)) = (f) (\9;u)) . (62) 
V J e\ He V JL / e V 1 e\ so 

in the right hand side we get no cross term since F is a symmetric function of 0's, see (47). 
Consider the cross-term in the left-hand side. Using the relations derived in the previous 
section we have 

L 

(f||0;u)) = -/^[P^ +1 ,F]P JJ+1 |0;u) +^[e Lil ,F]r u |0;u) (63) 

As usual we set the impurities to zero at the end in the right hand side. Note that at this 
order - since both terms are already multiplied by a g 2 factor - we can replace either of the 

kets |0;u) by ^|#;u}^ if we want to. Note also that only the last term in (63) does not 

have have nice cyclicity properties; this term can be interpreted as a correction to the state. 
More precisely, if we define 

|u) = (l-£ 2 r u e Lil )(|tf;u)) (64) 
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The eigenvalue reads 



L n , . ,„ N 



T(0;w)=a(0; ,)__ + __, where a(e;v) = [[ v _ 0a _. /2 , Q n (v) = [[(v-u j ) (59) 

a — 1 j — 1 
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then (62) becomes 

F(v)\u) = F(v)\u) (65) 

where 

L 

P = ( P )e ~ ^ D P W P]P ^ +1 and T = { F ) 6 ■ (66) 

The very exciting property we highlighted above is present again in (65). Namely, we have 
a set of states |u) which do not depend on the spectral parameter and still diagonalize a full 
family of operators J-(v) for arbitrary v. In this sense, we have integrability before and after 
the G-morphism. 



6.4 Two-loop Hamiltonian 



We found the first sign of evidence in favor of the G-morphism in (51). In these Bethe 
equations we identified the corrected momentum 

x(u + i/2) 



p(u) = - log . . 

1 x{u — i/2) 



Now we should look for the corrected energy 

2ig 



2ig 



x(u + i/2) x{u-i/2) 



(67) 



(68) 



More precisely we want to find some operator F such that P defined through (66) leads to 
the two loop Dilatation operator (2). We will look for this operator by fixing the eigenvalue 
T to be equal to the energy 



7 



E 



en, 



N 

E 



+ 



1\3 



+ c.c. 



Using (59) we can directly check that 
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? 2^ 2 cUogT + ^ 4 3>gT) 



u=-i/2 



(69) 



(70) 



Hence, following (66), the candidate operator that should correspond to the two loop Dilata- 

,14 



tion operator is 

L 

{2g 2 id u log 7 ) e + g A idl log f - 2g A ^[P JlJ+ i, id u log f]P jd+1 



(71) 



u=-i/2 



13 We would still get (69) without the 9-morphism Q , that is if we simply set the impurities to zero and 

not take any derivatives. Of course, in that case the rapidities would obey (37) with 8j = 0, that is they 
would be solutions to the one loop Bcthc equations which is not what we want. If we apply the 0-morphism 
we do get the two loop corrected Bethe equations (51) as desired. Indeed, as we can see directly from (72) 
below, we would not get the correct two loop Hamiltonian without the morphism contribution. 

14 At this loop order g 4 id^logf and g 4 idl(\ogf^ are indistinguishable. 
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u=-i/2 








+ 2g* 




-0 





Using the definition (57) and (52) we find indeed that these three contributions lead to 

(2g 2 id u \ogf) 1 ( 2g 2 -8g i 

gHdl logf 9 \ = \ - V \ £(I ~ P M+i) + 

+v 

-V J i=i 

Adding the three quantities we see that (71) is exactly equal to 

L L 

W = (2g 2 - 8<? 4 ) ^(1 - P M+1 ) + 2g 4 ^(1 - P i>i+2 ) (73) 

i=l i=l 

which is the Dilatation operator of M = 4 SYM up to two loops [8], see (2). Let us summarize 
the outcome of these last sections. We proved that 

The states (64), (31) are the precise two loop eigenvectors of the two loop Dilatation 
operator (2) of Af — 4 SYM. In particular, all the monstrous contact terms discussed in 
section 3 and appendix A are automatically incorporated in the algebraic approach (64) using 
the Q-morphism! 

This is one of the main results of this paper. Higher loops will be briefly touched upon in 
section 6.5. In section 7 we will apply these results to the computation of structure constants 
at one loop level. 



6.5 G— morphism and the All Loop Dispersion 

In this section we will briefly discuss the G-morphism from a more axiomatic point of view. 
The goal is to find out how we can constrain the B-morphism action at higher loops. We 
start by an ansatz of the form 



f(6)) = exp [J2J29 2l Mdi, ■ ■ -,d i+l ) f{6) 



,i=i i=i 



(74) 



where P21 is polynomial of degree 21. In what follows we assume that P21 is a homogeneous 
polynomial of degree 21. Note that in the exponent the polynomial only acts on / neighboring 
sites. We call this axiom the locality axiom. It is inspired by the structure of perturbation 
theory where the interaction range also increases with the perturbative loop order. 

Next we require the second axiom: 

(Um(O)h(O)) = (f sym (9)) (h(6)) . (75) 

which we call the morphism axiom. As explained before, this is a very important property 
which allows us to treat the Bethe roots as numbers, ignoring their dependence on the 
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impurities 9j. This seems like a very mild requirement but it turns out to be extremely 
powerful! We can study the constraint (75) in perturbation theory. We did it up to order g 8 
by considering polynomial functions of increasingly higher degree. The result of this brute 
force analysis is quite surprising. We found that the morphism axiom fixes completely 15 

0*). - 3 3 * 4 

(<?), = 359* (»?)," 35 (76) 

The action on 9? is not constrained and amounts to (re) defining the coupling constant so we 
set it to g 2 . We can ask the famous The On-Line Encyclopedia of Integer Sequences 16 



q2n 
3 , 



(77) 



which are very familiar in this context [8, 12]. We conjecture that (74) and (75) imply 
(77). We could not find a prove of this simple statement. If true, (77) immediately implies 
something quite remarkable: At all loop orders 



( 



log I 



u — 9i + \ \ x(u + 



» -L\og^ 1( (7* 



U — Vi — |/ X{U 



where x(u) are the Zhukowsky variables. We see that the one loop momentum (with impuri- 
ties) becomes promoted to the all loop momentum! We believe this is a very strong evidence 
in favor of this approach. We will continue our first principles study of the 0-morphism 
elsewhere [20]. 



7 G— morphism and One Loop Structure Constants 

In this section we will use the B-morphism to derive the one loop structure constants by 
promoting the tree level impure structure constants of section 5.1. Section 4 was also devoted 
to the study of one loop structure constants at one loop. The main differences between the 
approach in this section and the method used in section 4 are 

• First, here we will use the Algebraic Bethe ansatz approach together with the 0- 
morphism whereas in section 4 we followed the Coordinate Bethe ansatz. Of course, 

15 The morphism axiom is by no means sufficient to constrain (74) completely. It does not fix the form of 
the polynomials P21 at higher However, it does seem to fix the parameters of the polynomials that enter 
in the computations of ^0'- ^ . This is quite nice because it is exactly what we need to act on the dispersion 
relation! 

16 http://oeis.org/ or directly http://oeis.org/search?q=3%2C+10%2C+35 
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at the end of the day, both formalisms will give exactly the same result for physical 
quantities such as (the absolute value of) structure constants! 

However, in intermediate stages, scalar products in the coordinate and algebraic basis 
will differ since the kets are normalized differently. It is convenient to know the relative 
normalization to be able to compare not only the final results but also the intermediate 
ones. For example, at one loop we have (34). By checking a few examples like (32) 
at two loops, we can find the relative normalization between the two loop coordinate 
states - which we denoted by |p) in section 3 - and the algebraic states - which we 
denoted by |u) in section 6.3. We find 

|u> = fi |p> , H=U- ^Pj f[(e~ ip > - 1) HffaPi) (79) 

V 9 ' j=l i<j 

where (4) and (27). 

Another minor difference is that in the Coordinate Bethe ansatz (sections 3 and 4) we 
use mostly the momenta p a as variables. In the Algebraic Bethe ansatz (sections 5 and 
6) we use mostly the Bethe rapidities u a as our variables. The relation between the 
two is given by u a = u(p a ) where (10). Equivalently 

_ x(u + i/2) _ u+y/u*-4g2 _u_g_ 

6 ~ x(u-i/2)' X[U) ~ 2g ~g ' ' ' m 

We will also use x ± = x(u ± i/2) and x^ = x{u a ± i/2). For example, the above 
equation (79) can equivalently be written as 

N 



where 



N 



r " = E^XT' f(u,v)=(l + -!—) ( 1 + , 2 | i w 2 _i_ M ) ( 82 ) 



The set of momenta of the three operators 0%, O2 and O3 is, respectively k, p and q 
while the corresponding rapidities are denoted by u, v and w. 

The second important difference between the two approaches is in the method used. 
In the coordinate approach of section 4 we start with the scalar products for the one 
loop eigenstates and we guessed the higher loop generalizations. These conjectures 
were based on (1) playing with examples with a small number of excitations and find- 
ing patterns and (2) using the experience coming from higher loop corrections in the 
spectrum problem. Sometimes, one can guess not only the first quantum correction 
but even conjecture all loop results, see e.g. section 4.1. However, for the most compli- 
cated expectation value involved (6) the guessing strategy was manifestly insufficient 
to guess even the first quantum correction. With the algebraic Bethe ansatz approach 
we have full control over the two loop eigenvectors. As such, all one loop structure 
constants can be derived using the 0-morphism without any guessing. 

In sum, the two methods are nicely complementary. 
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7.1 Derivation of the Norm 



In (6) we identified all required scalar products and expectation values that we need to 
compute to obtain the complete one loop corrected structure constants in the SU(2) setup 
of [3]. We will now compute them using the algebraic Bethe ansatz approach. 

In this section we consider the simplest ones, the norms (1|1) etc. The norms were studied 
in the coordinate Bethe ansatz language in section 4.1. In the algebraic Bethe ansatz the 
corrected two loop eigenvectors are given by (64), |1) = |u) so what we want to compute is 

(u|u) = ((0;u|) fl (l-2p 2 r u EI L ,i)(|0;u)) e . (83) 

Note that for cyclic states the insertion of the operator EI will not depend of the position 
where it is inserted. This observation allows for the replacement H^i — > T u /L since we can 
average over where we insert. We will not need to use this in the derivation. 

The strategy is now to convert the contraction of two G-deformed states to the overall 
O-deformation of the scalar product 

((0;u|0;u)) . (84) 

Note that we do know the explicit expression for the quantity inside the brackets; it is given 
in (36). Hence evaluating (84) is straightforward (we will do it below). 

To bring (84) to a form closer to (83) we have to compute the cross-term. Using the 
identities derived in section 6.2, in particular (55), we easily find 

L 

(<0;u|][|0;u>) = ^.g 2 (0;u|(H ij+1 -^ L r u ) 2 |0;u) (85) 

= /(0;u|2r u + r 2 -2r u H liL |0;u) . (86) 

We see that the term with insertion of the Hamiltonian density is exactly the same as in 

(83) ! 17 Hence 

(u\u) = (l-2g 2 r u -g 2 r 2 u )((6;u\6;u)) , (87) 

V / 6 

which is exactly the kind of expression we were after. Furthermore, as mentioned above, 

(84) can be readily computed: the function (36) only depends on the impurities 9j through 
the symmetric combination appearing in (48). Hence, using (48) and (49) we immediately 
conclude that all we need to do is replace the dispersion factor, that is 



(<0;u|0;u) 



n ut - u m + i d 
11 , II >i h< N 



Uk - U m j,k<Ni dUj 
m?=k J 



L x{u k + i/2) lv^, u k -u m -i 
i x(u k -i/2) % u k -u m + i 



Taking into account the conversion factor (79) one can check that this is in perfect agreement 
with the conjecture for the norm of section 4.1! 

17 In (86) we can replace |0;u) by ^|#;u)^ if we want to. At this loop order they are indistinguishable 
since (86) already has an overall g 2 factor. 
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7.2 Derivation of the DVD Amplitude 

Having computed the norms we now move to the next line in (6), the contribution simple. 
As explained in section 4.2.1 this quantity can be simplified further so that all we need is 
the scalar product (19) between a Bethe state and a state obtained by gluing together two 
different vacuum descendents. 

For one loop eigenvectors with impurities (31) this quantity is given in (41). On the 
other hand, we are interested in the two loop eigenvectors with no impurities (64), that is 
(vac|v). To relate the two we follow the same strategy as in the previous section and start 
by applying the O-morphism to (41). We get 

^(vac|# (2) ;v)^ = (vac|(|0 (2) ;v)) = (vac|v) + g 2 T v (vac|H liL2 10; v) . 

The first term in the right hand side is exactly what we are after. As for the last one we can 
use (55) with j = L = L 2 to get 



vac 



v) = ((vac|# (2) ; v)) e - g 2 ^^ - id L , + T v ) (vac|^ 2 );v) (89) 



Note that in the right hand side of (41) there is no dependence on 9l 2 and thus we can drop 
the di 2 from (89). All we have to do now is compute the derivatives in the right hand side 
of (89). For that we use 

(et) = (1 - 9 2 rl)4 12 + 0(g 4 ) , g 2 id x (et) = ~9 2 T^ + 0(g 4 ) (90) 
where we use the short-hand notation 

x(a + V2) x i12 



As usual, the 0-morphism naturally gives rise to the Zhukowsky variables. The obtained 
result (vac|v) = (1 - # 2 r 2 ) £ (-1)^ f a ' a c^ 12 (1 + g 2 T a T a ) can be simplified further to 

aUc5=v 

(vac|v) = (l-£ 2 r 2 ) £ (-l)l*lp a ej", (92) 

aUa=v 

if we introduce a modified function 

f («, «) = ( l + T^t) ( l + 777o ■ i?L . ( 93 ) 



U-V)\ 1 ( U 2 + 1)( W 2 + 1) 

to absorb the partial energies r a and T & . The result we just derived coincides precisely with 
the conjecture in section 4.2.3 based on the coordinate Bethe ansatz approach! 18 



18 Of course, to compare both one needs to take (79) into account. 
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Using (92), (87), (51) we can write the final result (24) in the following succinct form 

(1-^/2) Yl (-i) |a| P*# a 

^yo.o yLiL 2 L 3 aUa=v , nA s 
C 123 - i = i > l y4 J 



which agrees perfectly with the previously conjectured expression in the coordinate Bethe 
ansatz approach of section 4. The case CI** is not harder to get, see 4.2.2. 



7.3 Two non-BPS Three Point Function C**3 and also C"' 

In this section we apply the technology developed above to the most complicated case of the 
structure constant when all three operators are non-BPS. The operator O3 can be replaced 
by the corresponding BPS operator without any loss of generality since its contribution 
factorizes and appears only in the simple contribution in (6). Hence, we will first consider 
this to be the case and recover the general case at the very end. 

Our starting point is again (64). We pick two closed states at two loops 

|1> = (l- J 2 r u H ll . 1 )(|9 ,1, ;u)) ()(i) (95) 

12) _ (1-/^,0(10^))^ 

Then we have to compute the most complicated contribution in (6) which is the last one, 
denoted as involved. We have 

involved = (1|4|2> - g 2 (l\ (H il2 , Ll2+1 + M LlA ) 6 3 (H Ll2jLl2+1 + M L2il ) |2> (96) 

where 

U=t4 iL 12 =t,i L 13 =N 3 L 23 =L 3 -N 3 

It is also convenient to define the one loop states as |1) = |0,u) and similar for |2). We 
shall not compute the second term in (96) since we will see soon that it neatly cancels with 
a contribution from the first one. Hence we focus on simplifying that first term 

(l|d 3 |2> = (((^1)^(63^^ (98) 

Again we have the G-morphism acting on the bra and ket states separately. The strategy is 
as before: we want to to bring it to the form of a differential operator in #'s acting on the 
tree level expression with impurities 

(0«;u|d 3 |0 (2) ;v). (99) 

This time it may seem a bit ambiguous since there are two different sets of #'s involved. 
However, as explained in figure 3, the scalar product (99) only depends on 9^ so that it 
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is rather natural to pick this set of 0's to act on. Moreover, the argument in that figure 
actually implies something slightly stronger, namely that the bottom part of the figure does 
not depend on 9^ and hence it only depends on 0^ C The bottom is nothing but 
3 \6^\ v) so we conclude that in (98) the first term is 19 



(<*»;u|) (<J,| S <V> 



0(2) 



u 



em 



O 3 |0 ( V> 



e(i) 



(100) 



By acting with the B-morphism on (99) we obtain the right hand side of this expression 
plus crossed terms that need to be massaged as usual. This is done in detail in appendix E. 
At the end of the day we find 



(i|d 3 |2) = (<0(Vl4|0 ( V>) +<7 2 (r v -r u ) 

V lew 



iV 



(i) 



r + r 

J- V ~ L u 



- (7 2 (r u + r v r u + r v )(i|03|2) 

+ g 2 (i\u Ll2A 6 3 - e Lljl 6 3 - o 3 H Ll2iLl2+1 - o 3 e L2il |2> . 



(0« u|Os|* ( V> 



(101) 



We see that the last line cancels precisely with the second term in (96)! Recall that these 
terms came from internal Feynman loop contributions. We see that they neatly cancel with 
particular contributions coming from the wave function corrections. 

More remarkable simplifications arise when we combine involved with the other contri- 
butions in (6) to obtain the physical quantity |C** 3 |- For O3 BPS we have 



I r^ m *° 1 

I *— ' 123 I 



' L1L2L 



L ^ x (involved/ V(l|l)(2|2) 



( Ls ) 

\N 3 J 



(102) 



Then dividing (101) by the norms of |1) and |2) given by (87) we get 

W;u|a 3 |^);v) NN 



involved/ v / (l|l)(2|2) 



'AM 



u^ 1 ); u)) ((0( 2 ); v|0( 2 ); v)) 



(103) 



0(2) 



+ 



g 2 (r v - r u ) 

V / (l|l)(2|2) 



v<£} - |(r v - r u ) (e^;n\d 3 \e^-v) 



As we show in Appendix E.l the last line is purely imaginary and thus can be absorbed into 
a physically irrelevant phase prefactor in the left-hand side! We conclude that the properly 

19 Very explicitly: 



(O 3 |0 (2) ;v>) ( 



= Os 



Ll g 2 



Os|v> + 



L 12 -l 



E to 



^(D +d 2 aW +d 2 0(1) 



6 3 |^;v)) e(i) = (6 3 |^;v)) e(2) 
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normalized structure constant is given by 



I r 1 ** 1 

l°123 I 



V L\L 2 L 3 




u|d 3 |* (2) ;v>] 


'em 






em 


r (0( 2 );v|^( 2 );v)^ 


'(9(2) 



(104) 



In sum, all ugly factors miraculously canceled out leaving us with the nicest possible expres- 
sion! 

From the discussion in section 2 we can easily generalize to the three BPS case which 
only differs by an additional factor corresponding to the third operator: 



I r* 19 " I 

l°123 I 



v ^X^((^ 1 );u|a 3 |^;v)) e(i) ^3(w)(l-^r w /2) 



(<0 (1) ;u|0W; u)) ew ((0W;v\eW; v)) m (3|3) 



(105) 



Finally, it is straightforward to compute all derivatives appearing in the several Q- 
morphisms since we know the impure results, see section 5.1. The final result for the one 
loop structure constants is presented in Appendix F. The expression (105) - or equivalently 
the more explicit expressions in the Appendix F - is the main practical outcome of the 
G-morphism. 



8 Conclusions 

Let us conclude with a small list of the main ideas/results together with a link to the 
corresponding equations. 

We constructed the eigenvectors of the two loop dilatation operator (2) both in Coordinate 
and Algebraic Bethe ansatz formalisms. Using the Coordinate Bethe ansatz this amounts 
to finding all contact terms, see (15) and (109). In the Algebraic Bethe Ansatz approach 
we used the G-morphism (43) to promote the one loop eigenvectors with extra parameters 
(impurities) (29) into exact two loop eigenvectors (64). In (17) we conjectured the form 
of the norm of the states at any loop order and made various tests of the conjectures by 
considering several states up to four loops. We also give a derivation using the algebraic 
approach at two loops in (87). 

Another important scalar product is the Double Vacuum Decay (DVD) amplitude (19), 
which we found up to two loops in (28) and (89). Using DVD as a building block we found 
the structure constants C^' (for O x BPS) in (24). 

Finally, the G-morphism techniques were used to derive the most general structure con- 
stants Cm (all are non-BPS) in (105) and (129). 

These expressions are basically ratios of simple determinants. As such they can be easily 
evaluated in a computer, even for operators with a large number of magnons which are the 
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relevant operators in the so called classical limit. This allows for nice comparisons with 
strong coupling which were discussed in the note [22] . 

It seems relatively straightforward to further correct the B-morphism and take into 
account in this way the all loop corrections to the dispersion relation [20]. In particular, 
we have full control over finite L states. On the other hand, it seems harder to incorporate 
the dressing factor. It would be very interesting to make contact with the boost operator 
approach of [9] where the dressing factor corrections are under control for large L states. To 
gain further insight into the underlying structure it is also very important to consider more 
general scalar products with fermions and derivatives. We hope that with some more data, 
a unifying structure will emerge and very general all loop expressions will become accessible. 
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A The four magnon contact term 

By directly diagonalizing the two loop dilatation operator (2) for large chains with four 
magnon excitations, we can derive the four particle contact term, see section 3. We find 
Ci„.(pi,p 2 ,P3,P4) = 2a /b where (y a = e tPa ) 

a = Zylvlvlvl - vlvlvlvl - vlvlvlvl + vlvlvlvl - vlvlvlvl + vlvlvlvl + vlvlvlvl 

%2 2 2 3 i 223,2 23 23,22 3 2 3 2 3 

2 2/s2/42/i + 2/22/32/42/1 + 2/22/32/42/1 - 2/22/32/42/1 + 2/22/32/42/1 - 2/22/ 3 2/42/i - 2/ 2 2/32/42/i 

3 3332, 2332, 3232 ^2232, 232, 2 32 

+ 2/22/32/42/1 - 2/22/32/42/1 + 2/22/32/42/1 + 2/22/32/42/1 - %22/ 3 2/42/i + 2/22/ 3 2/42/i + 2/ 2 2/32/42/i 

32, 222 22, 3322 2 3 2 2, 322, 222 a 3 2 2 2 

- 2/22/32/42/1 + 2/22/32/1 - 2/22/32/1 + 2/22/32/42/1 - 4?/22/32/42/l + J/22/32/42/1 + 2/22/42/1 - 4?/22/32/42/l 

1^2222 rr 222,222 22,3 22 t 2 2 2 i c 22 

+ I72/22/32/42/1 - 7J/22/32/42/1 + 2/32/42/1 - 2/22/42/1 + 2/22/32/42/1 - 7j/ 2 2/32/42/i + 62/22/32/42/1 

- 2/32/42/? - 2/22/32/? + 2/22/32/? + vlvlviVi ~ 2/22/32/42/? - 2/22/42/? + vlvlvivl ~ ^vlvlv^Vi 

+ 62/22/32/42/? - 2/32/42/? + 2/22/42/? - 2/22/32/42/? + 62/22/32/42/? - 72/22/32/42/? + 2/32/42/? + vlvlvlvi 

- 2/22/32/42/1 - 2/22/32/42/1 + 2/22/32/f2/i - 2/22/32/1 + 2/22/32/1 + vlvlvlvi - v-ivlvlvx - vlvlvi 

+ v\vlv\v\ ~ ^vlvlvlvi + 62/22/I2/42/1 - 2/32/42/1 + 2/22/42/1 - v\vzv\v\ + 62/22/32/42/1 - 72/22/32/12/1 

+ 2/32/42/1 + 2/22/1 + 2/22/32/1 - 42/22/32/1 + 2/32/1 - vlvlvm + Vivlvm + vlvm ~ vlvlvm 

+ 62/22/32/42/1 - 72/22/32/42/1 + vlvm ~ 42/22/42/1 + 2/22/32/42/1 - 72/22/32/42/1 + 172/22/32/42/1 - 42/ 3 2/42/i 

+ 2/42/1 - 2/1 + 2/22/32/4 - 2/22/32/4 - 2/22/32/4 + 2/22/32/4 - 2/2 + 2/22/3 - 2/3 - 2/12/12/4 + 2/22/32/4 + 2/22/4 

+ 2/22/32/4 - 42/22/32/4 + 2/32/4 - 2/4 + 3 

and 

j, 3333, 2233, 2323, 3223 r> 2 2 2 3 , 223, 2 23 

b = 2/22/32/42/1 + 2/22/32/42/1 + 2/22/32/42/1 + 2/22/32/42/1 - 22/22/32/42/1 + 2/22/32/42/1 + 2/22/32/42/1 

1 22 3, 3,2332,3232 r> 2 2 3 2 , 232,2 32,222 

+ 2/22/32/42/1 + 2/22/32/42/1 + 2/22/32/42/1 + 2/22/32/42/1 - 22/22/32/42/1 + 2/22/32/42/1 + 2/22/32/42/1 + 2/22/32/1 

1 3322 o2322, 322, 222 o3222,' 7 2222 /i 222 

+ 2/22/32/42/1 - 22/ 2 2/32/42/1 + 2/22/32/42/1 + 2/22/42/1 - 22/22/32/42/1 + 72/ 2 2/3 2/4 2/i - 42/22/32/42/1 

+ 2/32/42/? + 2/22/32/42/? - 42/22/32/42/? + 22/22/32/42/? + 2/22/32/? + vIvIvavI + vlvlvivl - ^vIvIvavI 

+ 22/22/32/42/? + 2/22/42/? + 22/22/32/42/? - 42/ 2 2/32/42/? + 2/32/42/? + vlvlvlvi + 2/22/32/42/1 + 2/22/32/1 

+ 2/22/32/42/1 + 2/22/32/42/1 - 42/22/32/42/1 + 22/22/32/42/1 + 2/22/42/1 + 22/22/32/42/1 - 42/22/32/42/1 + 2/32/42/1 

+ 2/22/1 + 2/22/32/1 - 22/22/32/1 + 2/32/1 + 2/22/32/42/1 + 2/22/42/1 + ^vlvlvm ~ ^vlvm + vlvm 

- 22/22/42/1 + 2/22/32/42/1 - 42/22/32/42/1 + 72/ 2 2/32/42/i - 2y 3 2/42/i + 2/42/1 + 2/22/32/4 + 2/22/32/4 + 2/22/3 

+ 2/22/32/4 + 2/22/4 + 2/22/32/4 - 22/22/32/4 + 2/32/4 + 1 



33 



B Simplifications for O x BPS 



In section 4.2.2 we explained that if 0\ is a BPS operator we can simplify considerably the 
most complicated contribution involved. Here we provide some more details. When 0\ is 
BPS, involved is given by 



i-g 2 



,L 12 +1 





I o I 




• • • 1 L 3 -N 3 \ 



l-2 2 (H Wl2+1 + H L2il )|2) 



where the dots stand for all possible intermediate states of length L 12 (over which we sum) 
We can drop the first two Hamiltonian insertions since they kill the bra. Then we use 



1 N1 1 






i ° i 




1 Ni-Ns 1 


| 


1 Lx \ 






• • • 1 L 3 -N 3 \ 


= < 


1 L!-N 3 \ 


1 L 3 -N 3 \ 



Finally we use the same argument as in section 4.2.1 to replace H £l2iLl2+1 +M.l 2 ,i by the full 
Hamiltonian which acts on the ket |2) to give the energy of the operator C?2. We conclude 

that for |1) = | | \ the involved contribution can be again written in terms of the DVD 

amplitude (19), see (23). 



C Code to generate contact terms 

Once we have guessed the form of the one loop scalar product (19) there is a nice trick to 
use it to derive all contact terms. Indeed 



A N (p) = ( | g || A I I p) = 2, ■ ■ ■ , N) (106) 



where i[)(ni, . . . ,un) is the wave function and L is the length of the Bethe state. Since all 
magnons are close together, it is given by 

^(1, 2, . . . , N) = (1 + <7 2 Cil( Pl , . . . , PN )) MP) exp ( J2 3P p A ( 107 ) 

p \i=i / 

where 

A(...,a,b,---) = Ua ~ Ub + * A(...,b,a,...), (108) 

Ua - U b - I 

and A(l, 2, . . . , M) = 1. Of course, u a = u(p a ) with u(jp) defined in (10). Hence, we have 

(M 
<n^:.Mn vs i = i (109) 

A.n{P) 

where An(p) is given by (28). In this way we obtained a closed expression for all two loop 
contact terms. It would be interesting to generalize this to higher loops, see next section and 
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In Mat hematic a, 20 

e[u_,L_:l]=((I+2 u)/(-I+2 u))~L (l+(32 I g~2 L u)/(l+4 u"2) ~2+0 [g] "4) ; 
f=(l+I/(#l-#2)) (l+g~2/((#l~2+l/4) (#2~2+l/4)))&; 

num[M_] :=Sum[A@@P Product [e [u[P [ [j] ] ] ,j] ,{j,l,M}] ,{P,Permutations@Range@M}] // 

{A [a___ , b_ , c_ , d___] : > (u [b] -u [c] +1 ) / (u [b] -u [c] -I ) A [a , c , b , d] / ; b>c , A@@Range [M] -> 

den[M_] : =Product [1/ (1/e [u [i] ] -l)Product [1/f [u[i] ,u[j]] ,{j,i+l,M}] ,{i,l,M}]* 

Sum [Product [-e [ j , M] Product [f [k , j] , {k , a [ [2] ] }] , { j , a [ [1] ] }] , 

{a, {Complement [u/@Range [M] ,#] ,#}&/@Subsets [u/@Range [M] ,{0,M}]}] ; 

Ct [M_] : =Ct [M] =l-num [M] /den [M] //SeriesCoef f icient [# , {g , , 2}] & 

For example, to get C^.(pi,p 2 ,P3) we run Ct [3] (and Simplify the result). 



D Details on A 

D.l Two Magnon Computation 

Here we derive (26) and its one loop quantum corrected analogue in detail. We have 

ri2— 1 



i P )= Yl ^K^ir^ut-.-ut •••>£' (no) 

i<7ii<7i 2 <L' m—i 

where the subscript L indicates the length of the state. The wave function is given by 

^(m,n 2 ) = ( e ipini+ip2n2 + e ip2ni+ipin2 S(p 2 , Pl )) x (111) 

(1 + S n ^ ni+l g 2 C i ^(p 2j pi) + S n2=L ^ ni=1 g 2 C^(p2,Pi)) . 

where we explicitly made manifest the last contact term contribution for the sake of clarity. 
Now 

^(p) = (HHh6H|p)= E ( 112 ) 

l<ni <n2<L 

Note that the sum goes only until L < L'. The only difference compared to (110) is the 
shortened length L < L' . Hence the last contact term in (111) is never triggered. Computing 
the sum (112) using (111) leads to (26) with f(p,k) replaced by f{p,k) in (27). Repeating 
the same exercise for a few cases with more magnons we find perfect agreement with (28). 



D.2 Fixing f from Impossibilities 

We have 



A L (p) = ( 



N 


1 


L 


L'-L \ 



20 Thc code becomes infinitely easier to read if converted into Standard Form. For that, copy the code into 
a Mathematica cell, select the cell, right click, and select "Convert To" and then "Standard Form". 
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where |p) is a Bethe state with rapidities p = {pi, . . . ,Pn} and length L'. To leading order 
this scalar product is given by (25). 

The scalar product Al({Pi, ■ ■ ■ ,Pn}) does not depend on L', directly. That is, it is sim- 
ply a function of the Bethe roots p and of the partial length L. 21 It yields the probability 
amplitude of finding the last V — L sites of the chain empty and the iV excitations homoge- 
neously distributed along the first L sites. For this probability amplitude to be non-zero we 
obviously need 

L > N . (114) 

Very interestingly, the leading order expression (25) - which was originally defined for L > N 
only - vanishes for L < N, see appendix C of [18]. This is a very nice fact which is of course 
not necessary a priori. It would be equally consistent - but not as nice - to have (25) for 
L > N only and simply define Al = for L < N. 

To leading order we have (25) and we note that 

SM = WJ) +(D{92) - (115) 

In this appendix we propose an alternative path towards obtaining the higher loop general- 
ization of (25). We make the following assumptions: 

• The scalar product continues to take the form (25) and we simply need to correct 
f(p,k) into f (p, k) = f(p,k) + 0(g 2 ). This holds provided the lengths involved are 
large enough so that wrapping interactions do not play a role. That is, for safety we 
assume that V 3> L in the usual asymptotic Bethe ansatz sense. 

• At any loop order we have (115), that is 



• We have 

A N -i(p 1 ,...,PN) = 0{g m ). (117) 

This is a generalization of Al{pi, ■ ■ ■ ,Pn) = for L < N which is a property of (25) to 
leading order as we discussed above. As mentioned above, this equality is some kind 
of a bonus that we did not have the right to ask for. As such, we will only require 
for a milder version of it to hold at higher loops, namely (117). We might also impose 
other similar conditions for other values of L < N but for the purpose of the current 
discussion we will only use (117). 22 

We will now show how these assumptions fix the next loop order completely. For that 
consider the N = 2 case. Using (116) we get rid of f(p2,Pi) in favor of f(pi,P2) in Al(pi,P2)- 

21 Of course, the momenta p do depend on L' since they are the solution to Bethe equations on a circle of 
length L' . So the scalar product does depend on L' in this sense. 

22 This condition somehow resembles an asymptotic condition, that is of absence of "wrapping" corrections. 
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Next we set L = 1 and impose Ai(pi,P2) = 0(9 4 ) following (117). This fixes f(pi,P2) to this 
loop order, see (26). We find 

I _|_ g*(Pl+P2) 

f(Pl,P2) = ip9 . g , T— (H8) 

which leads precisely to (27) which we obtained before by a completely different method! At 
the same time, A\(pi, . . . ,Pn) — 0{9 2 ) f° r N > 2. This confirms that at higher loops we can 
not ask for as much as for tree level. This is in agreement with the weaker condition (117). 

At the next few loop orders the conditions mentioned above seems to be valid. This 
statement is based on experimenting with configurations with small lengths and small number 
of magnons. For example, the condition A2(pi,P2,P3) = 0{g & ) seems to be indeed satisfied 
[20]. On the other hand, this condition by itself is not enough to fix f(p, k) completely. By 
playing with more impossible configurations, with other L's and N's, it might be possible to 
constraint this function much more, eventually to all loop orders. That remains to be seen. 



E Details of the derivation of the Structure Constant 



In this appendix we present the derivation of (101) in detail. As usual, the starting point is 
the action of the 0-morphism on an impure scalar product 

(<*«;u|03| e <V>X M = (^"'LO^^L + (<e (1, ;»l][°3|e<V>)jip) 

The first term in the right hand side is exactly what we are after, see (100), and the left- 
hand side is what we have already since we know (38). Hence we are left with the task of 
simplifying the cross-term. 

As before we can recast the derivatives into operators using (55). We use a shorthand 
notation V^f = d j a ) — d ' („) . We also omit 6's in bras and kets as always 6^ appears in the 

combination (9^; u| and similar for 6^ 



i=l 



Ui-\ 



(120) 



= ^ (1 \u|d32?W|v>+P«<u|6s2?«|v>+2?g<u|a32?g|v> 

i=l 
i=l 

where we use again that 03<9„(i) |v) = 3 <9 (2>|v) = for i > L\i- The first term can be 
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Figure 4: As an operator M.l 12 ,l 12 +i^^-l 12 ,Li2+i = ^£1, 1^3^2,1 = 0- The reason is that for 
the Hamiltonian density EI = I — P to give a non-zero result it must act on two sites where 
the spins do not point in the same direction. If that is the case for the bottom (that is when 
acting on O2) then it will not be the case when acting on the top (that is when acting on 
0\). Similarly, if it acts on two different spins in the top it will necessarily act on two equal 
spins in the bottom. These two cases are illustrated in the this figure. Hence we always get 
zero. 

simplified by converting the derivatives into operators 
£12— 1 £12— 1 £12— 1 

Vi(u\d 3 Vi\v) = (u|e M , +1 a 3 e M+1 |v> = ^ (u\w i>i+1 d 3 + o 3 e M+1 |v> (121) 

i=l i=l i=l 

Li L 2 

= (u| ^e vm a 3 + #3%+! - e il2iil2+1 a 3 - e Llil a 3 - a 3 e Ll2iLl2+1 - a 3 e L2il |v) . 

1=1 1=1 

We use that all the Hamiltonian densities acting on the state 3 give zero. Note that the first 
two terms are the complete one loop Hamiltonians acting on their eigenstates. Thus the first 
two terms are simply (u|(r u + T v )(9 3 |v). Let us simplify the remaining two terms in (120). 
Using (55) again we can rewrite them as 

(u|e Ll2 , Ll2+1 o 3 e il2 , Ll2+1 |v) + (u|(H w -r u )d 3 (H L2)1 -r v )|v) (122) 

= (u|r v r u a 3 - r v e Li a 3 - r u a 3 e i2 |v) , 

where we used that the terms with two EI in the l.h.s. vanish. The reason is that both 
Hamiltonian densities act with one leg on (9 3 . The Hamiltonian density is only nonzero 
when the sites on which it acts are different, however when contracted with O3 one of them 
receives J, whereas another one t and thus one of the operators always gives zero (see figure 
4). Thus the cross term finally gives 
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o*W) l = (u|r v r u o 3 - r v u LlA o 3 - r u o 3 e L2il |v) + 

9i 



(123) 



L 1 L2 

(u| ^H M+1 o 3 + ^o 3 e M+1 - e Ll2 , il2+1 d 3 -m Lu1 6 3 - d 3 u Ll2M2+1 - o 3 e L2>1 |v) . 

i=l i=l 

Combining all pieces together we get 

<i|3|2> — (<0«;u|d 3 |0 (2) ;v>) fl(i) = -/(v|(r u e Ll , 1 a 3 + r v a 3 e i2 , 1 )|v) (124) 

-/(u|r v r u d 3 - r v H Lljl d 3 - r u d 3 e L2il |v> 

V(u|(r u + r v )d 3 - u Ll2M2+1 d 3 - m Lul 6 3 - 6 3 m Ll2M2+1 - o 3 e L2il |v> 

which simplifies to 

(1|3|2) = ((^ 1 );u|a 3 |^ 2) ; v))^ (i) + .g 2 (r v - r u )(u|H Ll!l e> 3 - e> 3 H L2)1 |v} (125) 

- g 2 (u\(r u + r v r u + r v )d 3 - e Ll2il o 3 - u Lul d 3 - 6 3 u Ll2M2+1 - a 3 e L2il |v) . 

The second term in the r.h.s. can be written in terms of a total derivative: 

(u|H Zl)1 a 3 -a 3 H L2)1 |v) = (i(d w )+r u ( U |)a 3 |v) 

> L 12 L 12 + 1 ' 

+ (u\d 3 (i(d ffW -o)-r v |v) (126) 



1 L 12 v u 



(uia 3 |v) 



V;u\6 3 \eW;v)) n +.g 2 (r v -r l 



„--T)(l) _ r- , -n 
1 L 12 1 v t 1 u 



(u|e» 3 |v) 



Finally we get 
(1|3|2) = 

- £ 2 (u|(r u + r v r u + r v )d 3 - e Ll2il o 3 - e Lljl d 3 - o 3 e Ll2 , Ll2+1 - d 3 m L2A \v) 

Which is the expression used in section 7.3. 



E.l On an imaginary term 

In this section we discuss the contribution of the second line in (103): 



g 2 (T v - r u ) 



iV^l - \ (r v - T u ) (6U;u\6 3 \eW;v) , V 



Now, from (38) we see that 



(1) 

L\2 







L\2 



5(1) 



hi/2 







(127) 



il2+l 



(0;u|O 3 |0;v) ^ it' U 



Vn+i/2 

39 



x 



llm=llla=l v m -i/2 



x real 



(128) 



The derivative iDj^ can act on either of these three factors. After acting with the derivative 

we should set the impurities to zero since (127) already has an overal g 2 factor. If iD^ acts 
on the real function it gives a purely imaginary contribution which does not contribute to 
the absolute value of the structure constants. 23 If it acts on the second term it gives zero 
since the second term is a symmetric function of all the rapidities. Hence we are left with 
the actionof iT>l on the first factor of (128). That action leads to an imaginary part (which 
we drop as before) and to a real part. The real part cancels precisely with the V factors in 
(127)! Hence, the conclusion is simply that (127) does not contribute to the absolute value 
of the structure constants. 



F One Loop Structure Constants 

After all dust settles 24 we end up with the following explicit result for the one loop corrected 
structure constants: 

r»-^„ m \ - ^ L ^ L 3 I 1 ~ 9 2 (Tw + T 2 V - a„y)) Pjy 3 (u, v)„4jv 3 (w)£jy 3 (u, v) 4 
Cl23 (U ' V ' W) " B(w)B(w)B(v) +0{9h 

(129) 

Let us summarize all ingredients in this expression to render this appendix self-contained. 
The three operators are parametrized by 3 sets of Bethe roots: u = {%=i,...,jvi}, v = 
{vj=i,...,N 2 } an d w = {wj=i,...,Ar 3 }. Then, starting from the simplest contributions, we have 



1 w = E ,„2 , 1 ' Fuv = o ~2 7 1 ~9 Yl T^TI ' "uv^J^ ,,2! 1 _ S 



w — / , 2 1 1 ' uv — n / v 2 1 1 n / j 2 1 1 ' uv — / v 2 1 1 / v 2 1 1 



Then 



p * (u ' v) s 'm» f -. i+ i)n(«-^.) with q(2) s n 



and where the famous Zhuwkosky variables 



z - Uj 

3=1 



a; u + — - = - , z(u) = - - - + 0(g 3 ) . 
x{u) g g u 



23 Recall that the structure constants are defined up to a phase. Multiplying an operator O a by a phase 
does not affect its norm which is fixed by (O a Ol) but it changes the three point function by that phase. When 
bootstrapping higher point functions these phases cancel out anyway since summing over what is flowing 
amounts to inserting 1 = J2 a \®a)(O a \ in the higher point function. Hence the phase of the intermediate 
operator drops out as it should. 

24 That is, after the derivatives in the 0-morphism are evaluated 
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We now get to the interesting factors. We have 

e n f(M)n 



*4jv 3 (w) 



«Uo=w ada,ada adot 



x(a + i 



N3 



x(a — i/2) J 



ft ( x to - */ 2 > - A n f ( 

.7=1 v v 3 1 ' ' l<i<j<N 



(130) 



Wi,Wj) 



with 



f(w,z)=ll + 
Then we have the denominators 



1 + 



B(u) = 



\ 



n 



dot J-lilag'te+W 



S(u n , u m ) 

„. :m <J !<^<V 9«i ^ « - Z-/2) ^ I 



+y-\o g Uk - Ui - 1 



with similar expressions for B(u) and i5(w). The S-matrix 

_ u - W + 2 

£> (m, t>) = . 

v — it — i 

We are left with the last, most involved structure. We have 



S N3 (u, v) = P + <? 2 ((iV 3 + l)X? [+2] + (JV 3 - 1)X?[ +1 < +1 ] - 2a uv V^ 



where 



V EE 



2> 



h2l _ 



9T(«i) 


9T(v N2 ) 






97>i) 


8T(v N2 ) 






97>i) 


8T(v N2 ) 




<9mi 


07>i) 


8T(v N2 ) 


du Nl 




9T(«i) 


8T(v N2 ) 


<9mi 




0T(«i) 


8T(v N2 ) 








8T(v N2 ) 


<9mi 




07>i) 


8T{v N2 ) 




8u Nl 



q 2 (u 1 ) ... g7v 3 (Mi) gjv 3 +i(«i) 

q 2 (u Nl ) ... qN 3 (u Nl ) q Na+ i(u Nl ) 

q 2 {u 1 ) ... gjv 3 (Mi) giv a+ i +2 («i) 

q 2 (u Nl ) ... qN 3 (u Nl ) q Na+ i +2 (u Nl ) 

q 2 {ux) ... gjv 3 +i(Mi) gjv 3 +i+i(«i) 

g 2 («jv 1 ) ••• <?Ar 3 +i(-UjVi) gjv 3 +i+i(«ivi) 

g 2 («i) ••• ?iv 3 («i) gjv 3 +i+i(«i) 

q 2 {u Nl ) ... qN 3 (u Nl ) q Na+ i +1 (u Nl ) 



(131) 
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Q(u-i) fx(u-i/2)\ L Q(u + i) 
Q(u) + \x{u + i/2)J Q{u) 



• For symmetric configurations of Bethe roots a uv = and the above expressions simplify 
slightly. 

• If O3 is BPS the structure constant C*" reduces to C**3- This case is not considerably 
simpler than the general case. To obtain this case we replace T w — > 0, £>(w) — > J (^) 
and A Na (w) ^ 1 in (129). 

• If 0\ is BPS the structure constant C"* reduces to C°°°- This case is considerably 
simpler than the general case. The quantity Sn 3 does not appear in this case, only A 
type quantities arise. This case is discussed in the main text, see section 4.2.2. 

• The contribution (130) can be written as a determinant [23]. 
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